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Abstract.  An  expression  for  an  equilibration  time  scale  is  developed  that  characterizes  the  rate 
for  a group  of  test  particles  with  a Maxwellian  velocity  distribution  but  different  temperature 
from  that  of  Maxwellian  field  particles  to  relax  to  the  same  temperature  of  the  field  particles. 
The  expression,  which  can  be  used  for  any  value  of  the  Coulomb  logarithm,  is  based  on  onefold 
integral  expressions  for  Fokker-Plunck  velocity-space  friction  and  diffusion  coefficients  presented 
elsewhere  [Y.  Chang  and  C.  A.  Ordonez,  Phys.  Rev.  E 62,  8564  (2000)].  A comparison  with 
Spitzer’s  formula  shows  a noticeable  difference  when  the  Coulomb  logarithm  has  values  smaller 
than  10. 


INTRODUCTION 


An  equilibration  time  scale  is  presented  that  characterizes  the  rate  for  a group  of  test 
particles  with  a Maxwellian  velocity  distribution  but  different  temperature  from  that 
of  Maxwellian  field  particles  to  relax  to  the  same  temperature  of  the  field  particles. 
Spitzer  [1,2]  formulated  the  equilibration  time  scale  for  a two-component  plasma  based 
on  Chandrasekhar’s  formulas  for  the  velocity  space  friction  and  diffusion  coefficients 
[3,4].  Due  to  the  approximations  employed,  Spitzer’s  formula  for  the  equilibration  time 
scale  is  limited  to  plasmas  with  large  Coulomb  logarithm  values.  The  present  theory  is 
based  on  the  binary  collision  model,  but  is  not  restricted  in  applicability  by  the  value  of 
the  Coulomb  logarithm.  The  reader  is  referred  elsewhere  [5,6]  for  more  details  regarding 
the  applicability  of  the  theory  to  plasmas  confined  in  Penning  traps.  For  example,  the 
present  theory  can  be  used  to  predict  the  rate  at  which  an  antiproton  plasma  equilibrates 
with  a low-temperature  positron  plasma,  while  both  species  are  simultaneously  confined 
by  a nested  Penning  trap  [7], 
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THE  EQUILIBRATION  TIME  SCALE 

Consider  a single  test  particle  moving  in  a plasma  of  Maxwellian  field  particles.  The 
test  particle  will  exchange  its  energy  with  the  field  particles  due  to  collisions.  In  a single 
encounter  of  a test  particle  with  a field  particle,  the  exchange  of  energy,  A E,  is  given  by 

A E = (Aw2  + 2vAu||)  , (1) 

where  m is  the  mass  of  the  test  particle,  w is  its  speed,  Aw  is  the  magnitude  of  its  velocity 
change  due  to  the  collision,  and  Aw||  is  the  magnitude  of  its  velocity  change  along  its 
original  direction.  Averaging  Eq.  ( 1 ) over  a Maxwellian  velocity  distribution  for  the  field 
particles,  the  average  time  rate  of  change  of  the  test  particle  energy  is  written  as 

(A E)  = ((Aw2)  + 2w(At/||))  . (2) 

In  Eq.  (2),  (Aw||)  and  (Aw2)  are  the  usual  Fokker-Planck  velocity-space  friction  and 
diffusion  coefficients.  Let  us  suppose  that  a group  of  test  particles  have  a Maxwellian 
velocity  distribution  with  temperature  T,  which  is  different  from  that  of  the  field  particles 
of  temperature  Tf.  The  equilibration  time  scale  can  be  obtained  simply  by  averaging 
Eq.  (2)  over  a Maxwellian  velocity  distribution  for  the  test  particles.  The  rate  of  change 
of  the  test  particle  temperature  is  defined  through  the  expression, 

lnk^  = j(AE)fM(v,T)dv,  (3) 

where  n is  the  density  of  the  test  particle  species,  k is  Boltzmann’s  constant,  and/M  (v,T) 
is  the  Maxwellian  velocity  distribution  with  temperature  T, 

,"[v'T)  = n{^Yav(-iw'!‘)'  (4) 

Because  the  test  particle  velocity  distribution  is  isotropic,  the  integral  over  solid  angle 
can  be  carried  out  and  Eq.  (3)  can  be  written  as 

/”<*£>/«  (".JV*-  (M 

Substituting  Eq.  (4)  into  Eq.  (5)  and  carrying  out  a variable  change,  u — v/v th  with  the 
thermal  velocity  of  a field  particle  defined  as  vth  = yj 2kTj/mj , a non-dimensional  form 
of  Eq.  (5)  can  be  written  as 

dT  8^/2  roo 

1i  = iftkL  (A£>exP  (-Cu  )«  *‘,  (6) 

where  C = Tfm/(Tmf)  with  mf  the  mass  of  a field  particle.  The  next  step  for  the 
calculation  of  the  equilibration  time  scale  is  to  integrate  Eq.(6).  In  order  to  integrate 
Eq.  (6),  one  must  substitute  into  Eq.  (2)  expressions  for  friction  and  diffusion  coefficients. 
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SPITZER’S  APPROACH  TO  THE  EQUILIBRATION  TIME  SCALE 


Spitzer  employed  Chandrasekhar’s  expressions  for  the  friction  and  diffusion  coeffi- 
cients [2], 


4avth\G  (u) 

\Avi\)  = “ - 

T[) 

and 

2 = 2 (ai>ih)~  Aerf  (») 

n ,u 

where  a = 2 ft/m,  the  reduced  mass  is  /z  = mmj/{m  + mj),  A is  the  Coulomb  logarithm, 
the  time  scale  for  single  particle  interactions  is  r0  = (njvli,Trrlyl , ns  is  the  density  of 
the  field  particles,  the  interaction  radius  is  defined  as  r()  = ZZfC? / (87re0A:  (/ iTj/nif )), 
Z and  Zf  are  the  charge  state  of  a test  particle  and  field  particle,  e is  the  unit  charge,  f() 
is  the  permittivity  of  free  space,  erf  is  the  error  function,  and  Chandrasekhar’s  function 
is  defined  as  G (u)  = — (erf  (it)  /u).  Substituting  Eqs.  (7)  and  (8)  into  Eq.  (2),  the 
integral  in  Eq.  (6)  can  be  carried  out  as 

dT  4-m  (<»></,  )2  A C1/2(l  + C ~ 2«~' 1 ) , . 

di  'Ss/nn,  k (1+0'V2 


(7) 

(8) 


where  the  Coulomb  logarithm  has  been  approximated  as  being  constant.  Equation  (9) 
can  be  rearranged  into  the  standard  form  as 


dT  Tf-T 
dt  ~ T*!’il/er  ’ 


(10) 


where  r^J’ltzer  is  the  equilibration  time  scale  obtained  by  Spitzer  given  by 

Spitzer  _ 3m 7 ,ij  ( 47rf„  A2  fkT  kTf\'i,i 
cq  8\/2nnf\  \ZZfc- ) mf  j ' 


(11) 


Spitzer’s  equilibration  time  scale  applies  when  the  Coulomb  logarithm  is  large  both 
because  it  was  approximated  as  constant  and  because  an  approximation  was  employed 
in  the  process  of  obtaining  Chandrasekhar’s  expressions  for  the  friction  and  diffusion 
coefficients. 


NEW  APPROACH  TO  THE  EQUILIBRATION  TIME  SCALE 

Recently,  a variable  change  technique  (5,8-12]  has  been  developed,  which  has  been 
used  for  deriving  the  Fokker-Planck  coefficients.  Exact  onefold  integral  expressions  for 
Fokker-Planck  velocity-space  friction  and  diffusion  coefficients  have  been  obtained  from 
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the  new  technique.  These  expressions  make  it  possible  to  recalculate  the  equilibration 
time  scale  without  placing  a restriction  on  the  value  of  the  Coulomb  logarithm.  The 
expressions  for  the  friction  and  diffusion  coefficients  are  [5] 

avth  f /erf([7)  + erf(W)  exp(-(/2)  — exp(-VE2) 
t0u2  J \ us  V^us 


<Au|>  = - 


and 


(Aw2)  = 


{avthf  f erf((/)+erf(HQ^ 

tqu  J us 


(13) 


where  U = u+us,  W — u — us,  andu^  — A v/(av,h)  is  the  non-dimensional  variable  for 
velocity  change.  If  we  substitute  Eqs.  (12)  and  (13)  into  Eq.  (2),  Eq.  (6)  can  be  expressed 
as  the  following  twofold  integral: 


cLT  4 m(avth)2C3^‘2  f f / erf(CT)  + erf(W)  exp(-[/2)  — exp(— VC2)\  ududus 

dt  3y/nT0k  J J \ (2d-1  — l)-1  us  2-’v/7raw2  Jexp(Cu2)' 

(14) 


The  integral  in  Eq.  (14)  about  u can  be  further  reduced  to 


dT 

dt 


4m  {avti,)2  (1/2  (1  + C - 2a  *)  r / (u2s  \ 

3^r0k  (1  -E  C)3/2  J \ 1 + C/ 


us  1dus. 


Let  us  rewrite  Eq.  (15)  as 


dT  Tf  - T 
~dt  ~ 


(15) 


(16) 


where  the  new  equilibration  time  scale  is  obtained  from  Eq.  (16)  and  Eq.  (15)  as 


new 

eq 


3mm  fT 
16  H2 


^ J exp(- x)x  ldx 


-l 


(17) 


Here,  r = ( rifVavnr 2)  \ where  the  average  relative  speed  between  two  different  groups 
of  Maxwellian  particles  is  vav  = y8fcT'/(7r/r),  and  the  interaction  radius  is  defined  as  r = 
ZZje^K&irtokT').  The  non-dimensional  interaction  strength  is  defined  as  x — Hj{kT'), 
in  which  the  collision  strength  is  defined  as 


(18) 


where  A p — mAv  is  the  magnitude  of  the  momentum  transfer  for  a collision  event,  and 
the  reduced  temperature  is  defined  as 
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(19) 


r = 


+ 


No  approximations  are  used  to  arrive  at  Eq.  (17).  However  the  integral  in  Eq.  (17)  diverges 
due  to  the  long  range  Coulomb  interaction,  and  a cutoff  must  be  imposed  on  the  non- 
dimensional  collision  strength  x.  The  minimum  and  maximum  non-dimensional  velocity 
change  can  be  expressed  in  terms  of  the  average  relative  speed  and  the  Coulomb  logarithm 
as  [5] 


US.i 


U6,n 


— e~A 

»,k 

■ V(l  I? 

Vfh  ‘ 


(20) 

(21) 


For  the  two  temperature  system,  the  average  relative  speed  is  vav 
2 ((1  + C)/(K))1/2  vth>  and  the  integral  limits  for  Eq.  (17)  become 


^min  — ^ 

7T 

_ 4 

f'max  — * 

7T 


-2A 


Substituting  these  into  Eq.(l7),  we  obtain  the  equilibration  time  scale 

3 mm  f ( 47rf0  x 2 


new  _ 
'eq 


3 mnifT 


kTt 


a (kT  ^ 

8/i2r  (0,|e-2\  |)  “ 4v^F/(/r  (O.^e-w,^)  [zZje*)  [m  + mf 


(22) 

(23) 

3/2 

(24) 


where  the  zeroth-order  incomplete  gamma  function  F((),  .r,,,,,,,  :cmax)  is  defined  as 


e dx.  (25) 

min 

The  difference  between  Spitzer’s  approach  and  the  present  approach  is  characterized  by 
the  ratio 


--Spitzer 

eg 


Tnow 

'eq 


r(04e-2\i) 


2A 


(26) 


which  is  only  a function  of  the  Coulomb  logarithm.  With  this,  a comparison  of  Eq.  (24) 
and  Spitzer’s  result  Eq.  ( 1 1 ) is  shown  in  Fig.  I . The  new  result  Eq.  (24)  is  generally  larger 
than  the  traditional  one.  When  the  Coulomb  logarithm  is  larger  than  10,  the  difference  is 
less  than  5%.  However,  the  difference  can  be  substantial  when  the  Coulomb  logarithm  is 
less  than  ten. 
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